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ABSTRACT: 

A series of tests were conducted to measure the combustion performance of the Fastrac 
engine thrust chamber. During mainstage, the thrust chamber exhibited no large-amplitude chamber 
pressure oscillations that could be identified as low-frequency combustion instability or “chug”. 
However, during start-up and shutdown, the thrust chamber very briefly exhibited large-amplitude 
chamber pressure oscillations that were identified as chug. These instabilities during start-up and 
shutdown were regarded as benign due to their brevity. Linear models of the thrust chamber and the 
propellant feed systems were formulated for both the thrust chamber component tests and the flight 
engine tests. These linear models determined the frequency and decay rate of chamber pressure 
oscillations given the design and operating conditions of the thrust chamber and feed system. The 
frequency of chamber pressure oscillations determined from the model closely matched the frequency 
of low-amplitude, low-frequency chamber pressure oscillations exhibited in some of the later thrust 
chamber mainstage tests. The decay rate of the chamber pressure oscillations determined from the 
models indicated that these low-frequency oscillations were stable. Likewise, the decay rate, 
determined from the model of the flight engine tests indicated that the low-frequency chamber 
pressure oscillations would be stable. 


INTRODUCTION: 

NASA’s Marshall Space Flight Center (MSFC) was tasked with developing a 60,000 pound 
(267,670 N) thrust, pump-fed, LOX/RP-1 engine under the Advanced Space Transportation Program 
(ASTP). This government-led design has been designated the Fastrac engine, illustrated in figure-1. 

The X-34 b vehicle, illustrated in figure-2, was to use the Fastrac engine as the main 
propulsion system. The X-34 was to be a suborbital vehicle developed by the Orbital Sciences 
Corporation. The X-34 vehicle was to be launched from an L-1011 airliner. After launch, the X-34 
vehicle was to be able to climb to altitudes up to 250,000 feet (76,200 m) and reach speeds up to 
Mach 8, over a mission range of 500 miles (805 Km). The overall length, wingspan, and gross takeoff 
weight of the X-34 vehicle were to be 58.3 feet (17.8 m), 27.7 feet (8.4 m) and 45,000 pounds 
(200,752 N), respectively. 

This report summarizes the formulation of the linear model used in the study of the low- 
frequency combustion stability characteristics of the two Fastrac engine configurations. The first 
configuration was the Fastrac thrust chamber assembly (TCA) and its propellant feed system at 
MSFC's test stand 116. The second configuration was the Fastrac flight engine tested at NASA's 
Stennis Space Center (SSC). Also summarized is the Fastrac TCA design, the layout of the 
propellant feed system at MSFC's test stand 116, and the low-frequency combustion stability 
characteristics exhibited by the Fastrac TCA and propellant feed system during testing as was 
determined from high-frequency chamber pressure measurements. 

OBJECTIVE: 

One of the objectives of Fastrac TCA testing was to achieve a stable mainstage test that met 
a criteria for combustion stability as established by JANNAF standards 1 . This criteria was that the 
amplitude of the chamber pressure oscillations was to be 10% or less of the mean chamber pressure. 


“Approved for public release, distribution is unlimited. 
b The X-34 project was cancelled by NASA in March, 2001 . 
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Typically, in the context of high-frequency combustion instability, two stardards were to be 
met. First, the damp time of bomb-induced chamber pressure oscillations was to be 29 milliseconds or 
less. This damp time was determined from a combustion chamber acoustic frequency for the first- 
tangential (IT) mode of 1922 Hz. Second, the amplitude of the chamber pressure oscillations was to 
be 10% or less of the mean chamber pressure after the bomb-induced chamber pressure oscillations 
damped out. Both of these criteria were met during bomb testing whereby the high-frequency 
combustion instabilities were eliminated by effective acoustic cavity design 2 . 

In the context of intermediate-frequency combustion instability, a series of successful mainstage 
tests were conducted to measure the combustion performance of the Fastrac engine thrust chamber. 

The thrust chamber exhibited benign, yet marginally unstable combustion in that no damage to the thrust 
chamber assembly was ever sustained. A noticeble number of these tests had amplitudes of at least 
10%, and most of these tests had chamber pressure oscillations with a frequency on the order of 500 Hz. 
This frequency was too low to be identified as acoustic or high-frequency combustion instability and too 
high to be identified as low-frequency combustion instability. The source of the intermediate-frequency 
combustion instability or "buzz” was traced to the RP-1 venturi whose violently noisy cavitation caused 
resonance in the feedline downstream. Combustion was stabilized by increasing the throat diameter of 
the fuel venturi such that the cavitation would occur more quietly. Therefore, the marginally unstable 
combustion was related to the facilities and not to the thrust chamber itself 3 . 

Now, in the context of low-frequency combustion instability, the objective was to study the 
characteristics of the Fastrac TCA coupled with the propellant feed system to determine if low-frequency 
combustion instability or “chug” could occur. 



Figure-1: Fastrac Engine. 


Figure-2: X-34 Vehicle. 



TEST ARTICLE DESCRIPTION: 

As illustrated in figures 3 and 4, the components of the Fastrac TCA test article are the thrust 
chamber, nozzle, fuel manifold, LOX dome, injector faceplate, and acoustic cavities. For the TCA 
operating conditions, the combined LOX/RP-1 flow rate is about 197 Ibm/sec (89.56 Kg/sec) at a 
mixture ratio of 2.34. The chamber pressure is 650 psi (4.495 MPa). Details on the design of the 
acoustic cavities and the injector faceplate are discussed in related literature 2 . 

TEST FACILITY DESCRIPTION: 

Testing of the Fastrac engine TCA was performed at MSFC’s test stand 116. A typical test of 
the Fastrac engine TCA at test stand 116 is illustrated in figure-5. 
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Figure-3: Fastrac Thrust Chamber 
Assembly. 



SEAL GROOVES 


Figure-4: Fastrac Injector. 



Figure-5: A typical Fastrac TCA hot-fire test at MSFC test 
stand 116. 


The RP-1 feed system is presented in figure-6. The RP-1 feed system consisted of the RP-1 
tank, the RP-1 cavitating venturi, the RP-1 valve, the RP-1 injection manifold on the TCA, and 3- 
inch (7.62 cm) diameter RP-1 feed lines that connect all of these components together. The RP-1 
cavitating venturi had a throat diameter of 0.6121 inches (1.555 cm) and a discharge coefficient of 
0.96. The RP-1 valve was about 1.792 feet (0.546 m) long. Between the RP-1 cavitating venturi 
and the RP-1 valve was about 11.25 feet (3.429 m) of RP-1 feed line. Between the RP-1 valve and 
the RP-1 injection manifold on the TCA was about 8.959 feet (2.731 m) of RP-1 feed line. Not shown 
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in figure-6 were the “steerhorns” that started at about 6.376 feet (1.943 m) downstream of the RP-1 
valve. The steerhorns was a split or bifurcation in the RP-1 feed line. The bifurcated feed lines 
extended about 120° apart to send RP-1 to the two inlets on opposite sides of the RP-1 injection 
manifold. 



Figure-6: Test stand 116 RP-1 feed system. 



Figure-7: Test stand 116 LOX feed system. 
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The LOX feed system is presented in figure-7. The LOX feed system consisted of the LOX 
tank, the LOX cavitating venturi, the LOX valve, the LOX injection manifold on the TCA, and 3-inch 
(7.62 cm) diameter LOX feed lines that connect all of these components together. The LOX 
cavitating venturi had a throat diameter of 0.9 inches (2.286 cm) and a discharge coefficient of 0.96. 
The LOX valve was about 1.67 feet (0.509 m) long. Between the LOX cavitating venturi and the LOX 
valve was about 8.83 feet (2.691 m) of LOX feed line. Between the LOX valve and the LOX injection 
manifold on the TCA was about 3 feet (0.914 m) of LOX feed line. 


TEST PROCEDURE: 

A typical mainstage test was initiated by slightly opening the LOX valve to allow a small 
amount of LOX flow to chill the LOX feedline downstream and LOX injector prior to ignition. During 
pre-chill, this small amount of LOX flow vaporized upon injection into the chamber and was 
discharged out through the nozzle. Ignition was performed by injecting TEA/TEB into the chamber 
while maintaining this small amount of LOX flow. After ignition, fuel-rich “level-1” combustion was 
initiated by completely opening the RP-1 valve to allow the full amount of RP-1 flow to be injected 
into the chamber. Simultaneously, the TEA/TEB flow was terminated and the small amount of LOX 
flow was maintained. After level-1 combustion was briefly established, oxidizer-rich “level-2” 
combustion was initiated by completely opening the LOX valve to allow the full amount of LOX flow to 
be injected into the chamber. Level-2 combustion constituted mainstage operation of the Fastrac 
TCA. 

Chamber pressures for level-1 and level-2 combustion were 450 psi (3.06 MPa) and 650 psi 
(4.42 MPa), respectively. 

Those feed systems inactive during pre-chill, ignition, level-1, and level-2 operations were 
purge with nitrogen gas to prevent backflow into the injector. 

Table-1: Fastrac TCA Mainstage Test Results 


Test # 

Type 

HHSHi 


■ 

IHKSnHH 

1 

ignition 

n/a 

n/a 

n/a 

n/a 

2 

ignition 

n/a 

n/a 

n/a 

n/a 

3 

level-1 

5 

410 

450/3.06 

11.8 

4 

level-2 

7 

513 

644/4.38 

5.1 

5 

level-2 

25 

508 

625/4.25 

11.1 

6 

not successful 

n/a 

n/a 

n/a 

n/a 

7 

level-2 

25 

522 

643/4.37 

8.1 

8 

level-2 

53 

522 

650/4.42 

10.2 

16 

level-2 

8 

527 

641/4.36 

9.2 

17 

level-2 

27 

527 

631/4.29 

10.5 

19 

level-2 

130 

522 

640/4.35 

11.9 

20 

level-2 

17 

508 

645/4.39 

8.5 

21 

level-2 

17 

503 

641/4.36 

8.9 

22 

level-2 

17 

493 

635/4.32 

12.9 

23 

level-2 

150 

498* 

656/4.46 

3.0 

24 

level-2 

150 

515* 

664/4.52 

3.9 

25 

level-2 

6 

142 

649/4.41 

6.8 

26 

level-2 

150 

★ 

664/4.52 

2.9 


SUMMARY OF MAINSTAGE TESTS: 

The results of the mainstage tests are presented in table-1. The amplitudes presented were 
determined by visual averaging and are not RMS values, which will be lower. Out of the 15 successful 
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mainstage tests, 10 of these tests had amplitudes of at least 8%, and 6 of these tests had 
amplitudes of at least 10%. These noticebly large-amplitudes in chamber pressure were mostly 
associated with the 500 Hz intermediate-frequency combustion instability or buzz. 

Tests 1 and 2 were ignition tests. Test-3 was a level-1 test that exhibited large amplitude 
chamber pressure oscillations at 410 Hz. Test-6 was unsuccessful. Tests 9-15, and 18 were not 
included in table-1 since these were bomb-tests 2 . Tests 16, 17, and 19 were level-2 tests conducted 
to obtain acoustic cavity temperatures. Tests 23 and 24 were level-2 tests that not only had the 500 
Hz oscillation, but had oscillations at 60 Hz and its harmonics indicated by the asterisk in the 
frequency column. Test-25 was the first level-2 test where the 500 Hz oscillation was eliminated, 
leaving behind what could be interpreted as low-amplitude chug. Test-26 was the second level-2 test 
where the 500 Hz oscillation was eliminated leaving only low-amplitude oscillations at 60 Hz and its 
harmonics indicated by the asterisk in the frequency column. Tests 27-31 also were not included in 
table-1 since these were additional bomb tests. 

CHARACTERIZATION OF CHAMBER PRESSURE OSCILLATIONS: 

Presented in table-2 are the natural frequencies of the combustion chamber acoustic modes. 
These frequencies were based on a chamber diameter of 13.28 inches (33.73 cm) and a chamber 
speed of sound of 3628.8 ft/sec (1106.1 m/sec). The frequencies of chamber pressure oscillations of 
interest were too low to be high-frequency. Also, the frequencies of chamber pressure oscillations of 
interest were lower than the 500 Hz buzz that was traced to the fuel venturi. Therefore, the frequency 
of the chamber pressure oscillations of interest were too low to be considered intermediate- 
frequency. Therefore, the chamber pressure oscillations of interest were low-frequency. 

Table-2: Fastrac TCA Chamber Acoustic Frequencies 


Tangential modes Radial modes Mixed modes 


IT 

1922 Hz 

1 R 

4000 Hz 

IT-1 R 

5565 Hz 

2T 

3188 Hz 

2R 

7322 Hz 

1T-2R 

8900 Hz 

3T 

4385 Hz 

3R 

10,618 Hz 

2T-1R 

7000 Hz 

4T 

5550 Hz 






REVIEW OF MODEL IMPLEMENTATION FOR THE PRESSURE-FED FASTRAC TCA: 

To perform qualitative predictions of low-frequency combustion instability, a linear model was 
formulated. The formulation of this model was based on the model of Wenzel and Szuch 4 with two 
very important improvements. The first improvement was the modeling of each of the propellant feed 
systems extending from the injectors up to the venturis. The second improvement was the 
assumption of a non-zero damping coefficient. Therefore, not only was it possible to generate a 
neutral stability boundary for a damping coefficient of zero, it was also possible to determine the 
frequency and damp time for chamber pressure oscillations given the steady-state chamber pressure 
and steady-state injector pressure drops. The details on the formulation are lengthly and will not be 
discussed here, but in the appendix. 

The model was implemented by first specifying the parameters of the Fastrac TCA 
combustion chamber, the oxidizer feed system, and the fuel feed system, all at test stand 116 during 
test-8. 


For the combustion chamber, the volume from the injector faceplate down to the throat was 
determined from the formula of characteristic length 5 : 


L * - “ - L cy £ c + “( £ c + + X ) • 

( 1 ) 
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Where A T is the throat cross-sectional area, Lev is the combustion chamber length from the injector 
face to the nozzle inlet, Leo is the nozzle length from the inlet to the throat, and e c is the chamber 
contraction ratio. The from figure-3, Fastrac TCA combustion chamber was 13.28 inches (33.73 cm) 
diameter at the injector face, and 8.25 inches (20.96 cm) diameter at the throat. This gave a 
contraction ratio and a throat cross-sectional area of 2.59 and 53.45616 in 2 (344.88 cm 2 ), 
respectively. Additionally, the distance from the injector face to the nozzle inlet was 2.7 inches (6.86 
cm), and the distance from the nozzle inlet to the throat was 15.3 inches (38.86 cm). Substituting 
these values into equation-1 gave a characteristic length of 33.52018 inches (85.14 cm). Multiplying 
the characteristic length and the throat cross-sectional area gave a chamber volume of V c = 1.03696 
ft 3 (0.02936 m 3 ). 

For a chamber pressure of 650 psi (4.422 MPa) and a mixture ratio of 2.34, a thermochemical 
equilibrium calculation determined that the temperature and density of the combustion products was 
Tc = 6397°R (3550°K) and pc = 0.212 lbm/ft 3 (3.403 Kg/m 3 ), respectively. Since the total flow rate was 
197 Ibm/sec (89.55 kg/sec) and the combustion chamber flush/fill time was defined as: 



( 2 ) 

The combustion chamber flush/fill time was Tc = 1.11592 x 10 3 seconds. Also, at a chamber pressure 
of 650 psi and a mixture ratio of 2.34, a thermochemical equilibrium calculation determined that the 
derivative of the chamber gas temperature with respect to the mixture ratio was -228.29°R 
(-126.67°K). The test-8 chamber pressure was 658 psi (4.476 MPa). 

The cross-sectional areas of the opened LOX valve and the LOX dome inlet were 7.06858 
in 2 (45.60 cm 2 ) and 4.53144 in 2 (29.24 cm 2 ), respectively. These were based on local feedline 
diameters of 3.0 inches (7.62 cm) and 2.402 in (6.10 cm), respectively. However, the total cross- 
sectional area of all of the LOX orifices was 3.04797 in 2 (19.66 cm 2 ). The lengths of the LOX valve, 
the LOX dome inlet and the faceplate thickness for the LOX orifices were 20 in (50.80 cm), 6 in 
(15.24 cm) (assumed), and 0.24250 in (5.08 cm). 

The volumes of the LOX feed system segments were 748.98710 in 3 (12,273.70 cm 3 ), 
163.13171 in 3 (2673.25 cm 3 ), and 287.41303 in 3 (4709.86 cm 3 ) for the LOX feedline upstream of the 
valve, the LOX feedline upstream of the injector inlet, and the LOX dome, respectively. 

The LOX feed system pressures were 1088.1 psi (7.402 MPa), 936.3 psi (6.369 MPa), and 
784.5 psi (5.337 MPa), upstream of the LOX valve, upstream of the LOX dome inlet, and in the LOX 
dome, respectively. All of the LOX feed system pressures were obtained from test-8. 

The LOX flow rate was 138 Ibm/sec (62.73 kg/sec) and the time lag required for 50% of the 
LOX droplet mass to evaporate was previously determined to be 3.00775 x 10^ seconds 6 . 

The cross-sectional areas of the opened RP-1 valve and the RP-1 injector manifold inlet were 
5.49467 ft 2 and 5.49292 ft 2 , respectively. These were based on feedline diameter of 2.645 in (6.72 
cm) upstream of the “steerhorn” bifurcation and 1.870 in (4.75 cm) downstream of the steerhorn 
bifurcation, respectively. However, the total cross-sectional area of all of the RP-1 orifices was 
1.35362 in 2 (8.73 cm 2 ). The lengths of the RP-1 valve, the RP-1 injector manifold inlet and the 
faceplate thickness for the RP-1 orifices were 21.50400 in (54.62 cm), 6 in (15.24 cm) (assumed), 
and 0.24250 in (0.62 cm). 

The volumes of the RP-1 feed system segments were 954.25877 in 3 (15,637.50 cm 3 ), 
590.66630 in 3 (9679.29 cm 3 ), and 89.70048 in 3 (1469.93 cm 3 ) for the RP-1 feedline upstream of the 
valve, the RP-1 feedline upstream of the injector inlet, and the RP-1 injector manifold, respectively. 
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Note that the cross-sectional area of the RP-1 injector manifold inlet and the volume of the 
RP-1 feedlines upstream of the RP-1 injector manifold reflect the actual bifurcated “steerhorns” 
configuration of feedlines not illustrated in figure-6. 

The RP-1 feed system pressures were 832.30 psi (5.662 MPa), 810.25 psi (5.512 MPa), and 
788.20 psi (5.362 MPa), upstream of the RP-1 valve, upstream of the RP-1 injector inlet, and in the 
RP-1 injector manifold, respectively. All of the RP-1 feed system pressures were obtained from test-8. 

The RP-1 flow rate was 59 Ibm/sec (26.82 kg/sec) and the time lag required for 50% of the 
LOX droplet mass to evaporate was previously determined to be 3.58088 x 10' 2 seconds 6 . 

The low-frequency combustion stability map for the Fastrac TCA with the pressure-fed 
propellant feed system is presented in figure-8. The test-8 RP-1 injector pressure drop was 0.19787. 
The test-8 LOX injector pressure drop was 0.19225. 

To generate a curve on the low-frequency combustion stability map, first a frequency 
bandwidth was selected. In figure-8, a typical bandwidth coresponding to zero growth rate was from 
67 Hz to 103 Hz. Secondly, a pressure oscillation growth rate was selected. A positive growth rate 
indicated unstable oscillations in pressure. A negative growth rate indicated stable oscillations in 
pressure. A zero growth rate indicated neutrally stable oscillations in pressure. The growth rate was 
arbitrarily varied from -75/sec to 75/sec. 

According to equation-129, solving for the RP-1 injector pressure drop required solving a 
quadratic equation. Since equation-129 was solved by implementing the quadratic formula in 
equation-130, both “positive” and “negative” solutions for the RP-1 injector pressure drop were 
searched for the actual solution. Then, equation-128 was used to determined the LOX injector 
pressure drop from the RP-1 injector pressure drop. By trial and error, a growth rate was selected until 
the curve passed through the injector operating point. Therefore, the growth rate of the amplitude 
and the frequency of pressure oscillations were determined in this manner. Therefore, in figure-8, the 
frequency and growth rate of pressure oscillations in test-8 of the Fastrac TCA, tested at test stand 
116 were 64.93 Hz and -68.5/sec (stable), respectively. 

Interpreting the frequency was straight forward. However, interpreting the growth rate was not 
as trivial. If one assumed that the damp time was the amount of time required by the pressure 
oscillation to damp down to 50% of its original amplitude, then: 

(0.69315) 

^DAMP _ ^ 

(3) 


Therefore, with a growth rate of -68.5/seconds, a pressure oscillation of 64.93 Hz would damp to 
50% of its initial amplitude in 10.12 milliseconds. In other words, a mainstage pressure oscillation with 
a frequency of 64.93 Hz and an initial amplitude of 20% of the mean chamber pressure would damp 
to an amplitude of 10% of the mean chamber pressure after 10.12 milliseconds. This damp time was 
extremely short if one considered that for the maximum damp time 1 , given by: 

(l250 msec-Hz 1/2 ) 
t 

l "MAX damp 

(4) 


4co 


gave 155.13 milliseconds. 

Also, in figure-8, it was apparent that the first fictitious operating point APf/Pc = 0.200 and 
APo/Pc = 0.126 was an intersection of several curves. This was a situation where a single 
combination of APf/Pc and APo/Pc gave several combinations of frequencies and growth rates. 
Recall that a single frequency and growth rate represents a damped (for negative growth rate), 
sinusoidal function, which was a solution of a set of ordinary differential equations. Since there were 
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no initial conditions given, all possible solutions will stand ambiguously. The only certainty for such a 
point is that since the growth rates for all three curves were negative, then the pressure oscillations 
would be stable regardless of the frequency. 

Additionally, in figure-8, it was apparent that the second fictitious operating point AP F /Pc = 
4.500 and AP 0 /Pc = 0.060 was an intersection of a couple of curves. This was a situation where a 
single combination of AP F /Pc and AP 0 /Pc gave two combinations of frequencies and growth rates. 
Again, recall that a single frequency and growth rate represents a damped (for negative growth rate), 
sinusoidal function, which was a solution of a set of ordinary differential equations. Since there were 
no initial conditions given, all possible solutions will stand ambiguously. Unlike the first fictitious 
operating point, the stability would be less certain since one curve had a positive growth rate and the 
other curve had a negative growth rate. Although this operating point was more likely to be unstable, 
it would not be a realistic operating point due to the extremely high RP-1 injector pressure drop. 
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Figure-8: low-frequency combustion stability map for test-8 of the 
pressure-fed Fastrac engine TCA tested at MSFC test stand 116. 


REVIEW OF MODEL IMPLEMENTATION FOR THE PUMP-FED FASTRAC TCA: 

The model was again implemented by first specifying the parameters of the Fastrac TCA 
combustion chamber, the pump-fed oxidizer feed system, and the pump-fed fuel feed system, all 
integrated together as a flight engine, tested at SSC. 

For the combustion chamber, the geometry, the chamber pressure, the mixture ratio, and the 
total flow rate are all unchanged in the conversion from a pressure-fed to a pump-fed propellant feed 
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system. The only changes are in the LOX and RP-1 feed systems. These changes in the LOX and 
RP-1 feed systems were reflected in the line length and line diameter measurements in figures 9 and 
10, repectively. Unfortunately, the venturis were not replaced by a pumps as a flow rate-setting 
device. 




Figure-9: LOX feedline segments on the upper Figure-10: RP-1 feedline segments on the upper 
portion of the Fastrac flight engine (front view). portion of the Fastrac flight engine (left view). 

Referring to figure-9, the lengths and diameters along the LOX feedline are: 

L x =■ 6.50 in (16.51 cm), D 3 = 2.645 in (6.72 cm) 

L 2 = 4.33 in (11.00 cm), D 2 = 2.645 in (6.72 cm) 

L 3 = 31.9 in (81.03 cm), D 3 = 2.402 in (6.10 cm) 

Referring to figure-10, the lengths and diameters along the RP-1 feedline are: 

L x = 19.51 in (49.56 cm), D 3 = 2.635 in (6.69 cm) 

L 2 = 4.33 in (11.00 cm), D 2 = 2.645 in (6.62 cm) 

L 3 = 27.46 in (69.75 cm), D 3 = 1.870 in (4.75 cm) 

The cross-sectional areas of the opened LOX valve and the LOX dome inlet were 
unchanged at 5.49467 in 2 (35.45 cm 2 ) and 4.53144 in 2 (29.24 cm 2 ), respectively. These were based 
on local feedline diameters of 2.645 inches and 2.402 in, respectively. However, the total cross- 
sectional area of all of the LOX orifices was 3.04797 in 2 (19.66 cm 2 ). The lengths of the LOX valve, 
the LOX dome inlet and the faceplate thickness for the LOX orifices were 4.33 in (11.00 cm), 6 in 
(15.24 cm) (assumed), and 0.24250 in (0.62 cm). 

The volumes of the LOX feed system segments were 35.71532 in 3 (585.27 cm 3 ), 144.55282 
in 3 (2368.80 cm 3 ), and 287.41303 in 3 (4709.86 cm 3 ) for the LOX feedline upstream of the valve, the 
LOX feedline upstream of the injector inlet, and the LOX dome, respectively. While the LOX dome 
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volume was unchanged, the volumes of the LOX feedlines upstream of the valve and upstream of 
the LOX injector inlet were reduced to reflect the shortening of the LOX feedlines. 

The LOX feed system pressures were unchanged at 1088.1 psi (7.402 MPa), 936.3 psi 
(6.369 MPa), and 784.5 psi (5.337 MPa), upstream of the LOX valve, upstream of the LOX dome 
inlet, and in the LOX dome, respectively. All of the LOX feed system pressures were obtained from 
test-8. 


The LOX flow rate was unchanged at 138 Ibm/sec (62.73 kg/sec) and the time lag required 
for 50% of the LOX droplet mass to evaporate was also unchanged at 3.00775 x 10' 2 seconds. 

The cross-sectional areas of the opened RP-1 valve and the RP-1 injector manifold inlet were 
unchanged at 5.49467 ft 2 (0.51047 m 2 ) and 5.49292 ft 2 (0.51031 m 2 ), respectively. These were 
based on feedline diameter of 2.645 in upstream of the “steerhorn" bifurcation and 1.870 in 
downstream of the steerhorn bifurcation, respectively. However, the total cross-sectional area of all of 
the RP-1 orifices was 1.35362 in 2 . The lengths of the RP-1 valve, the RP-1 injector manifold inlet and 
the faceplate thickness for the RP-1 orifices were 4.33 in (11.00 cm), 5.36 in (13.61 cm), and 
0.24250 in (0.62 cm). 

The volumes of the RP-1 feed system segments were changed to 106.39186 in 3 (1743.45 
cm 3 ), 178.14401 in 3 (2919.26 cm 3 ), and 89.70048 in 3 (1469.93 cm 3 ) for the RP-1 feedline upstream of 
the valve, the RP-1 feedline upstream of the injector inlet, and the RP-1 injector manifold, 
respectively. 
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Figure-11: low-frequency combustion stability map for the pump-fed Fastrac engine TCA tested at SSC. 
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Figure-12: low-frequency combustion stability map for the pressure-fed Fastrac engine TCA 
tested at MSFC and the pump-fed Fastrac engine TCA tested at SSC. 


As in the case of the pressure-fed TCA, note that the cross-sectional area of the RP-1 
injector manifold inlet and the volume of the RP-1 feedlines upstream of the RP-1 injector manifold 
reflect the actual bifurcated “steerhorns” configuration of feedlines not illustrated in figure-6. 

The RP-1 feed system pressures were 832.30 psi (5.662 MPa), 810.25 psi (5.512 MPa), and 
788.20 psi (5.362 MPa), upstream of the RP-1 valve, upstream of the RP-1 dome inlet, and in the 
RP-1 injector manifold, respectively. All of the RP-1 feed system pressures were obtained from test-8. 

The RP-1 flow rate was unchanged at 59 Ibm/sec (26.82 kg/sec) and the time lag required for 
50% of the LOX droplet mass to evaporate was 3.58088 x 10' 2 seconds. 

The low-frequency combustion stability map for the Fastrac TCA with the pump-fed propellant 
feed system is presented in figure-11. The LOX and RP-1 injector pressure drops for the SSC flight 
engine were assumed to be the same as the test-8 LOX and RP-1 injector pressure drops, which was 
0.19225 and 0.19787, respectively. A typical frequency bandwidth corresponding the zero growth 
rate was from 113.8 Hz to 137.5 Hz. Frequency decreased as APf/Pc increased. This was the exact 
opposite of figure-8, where frequency increased as AP F /Pc increased. As in figure-8, a pressure 
oscillation growth rate was selected. A positive growth rate indicated unstable oscillations in 
pressure. A negative growth rate indicated stable oscillations in pressure. A zero growth rate 
indicated neutrally stable oscillations in pressure. The growth rate was arbitrarily varied from -75/sec 
to 75/sec. Along the neutral stability curve, AP 0 /Pc seemed to be constant with AP F /Pc until AP F /P C 
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was about 100%, Below a AP F /Pc of 100%, AP 0 /Pc decreased rapidly with decreasing AP 0 /Pc- 
Likewise, the curves corresponding to positive growth rate and negative growth rate followed the 
trends of the neutral stability curve as they ran below and above it, respectively. Despite all trial-n- 
error attempts to choose a negative growth rate such that would produce a curve that would run 
through the injector operating point, none could be found. The best curve that could be produced 
corresponded to a growth rate of -128/sec. Therefore, the negative growth rate and the frequency of 
pressure oscillations corresponding to the injector operating point could not be determined. All that 
could be said about this point was that the pressure oscillations were likely to be qualitatively stable. 

In figure-12, the neutral stability curves for the pressure-fed and the pump-fed Fastrac TCA 
are compared. The curve for the pressure-fed TCA was obtained from the “negative” solution of 
equation-130. The curve for the pump-fed TCA was obtained from the “positive” solution of equation- 
130. For the curve corresponding to the pressure-fed TCA, the frequency increased as AP F /Pc 
increased. However, for the curve corresponding to the pump-fed TCA, the frequency decreased with 
increasing AP F /P C . By comparing the behavior of the neutral stability curves, one can conclude that 
for extremely high valves of AP F /P C , shortening the propellant feedlines to convert from a pressure- 
fed system to a pump-fed system will make no difference in the low-frequency combustion stability of 
the TCA. However, one can conclude that for realistically low values of AP F /P C , shortening the 
propellant feedlines to convert from a pressure-fed system to a pump-fed system will make a dramatic 
increase in the low-frequency combustion stability of the TCA. 



time( seconds) 


Figure-13: test-23 chamber pressure. 


REVIEW OF TEST DATA 

In figure-13, the test-23 chamber pressure is presented as a function of time. This chamber 
pressure vs. time data is also typical of tests 24, 25, and 26. With the exception of test-25, these 
tests were 150 seconds in duration. In these tests the chamber pressure data was divided into four 
phases. The first phase was start-up. The second phase was buzz-dominated mainstage. The third 
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phase is buzz-minimized or buzz-free mainstage. The fourth stage is shutdown. Tests 23 and 24 had 
all four phases. However, tests 25 and 26 had only three phases since there was no buzz-dominated 
mainstage. 

In figure-14, the test-23 chamber pressure during the start-up phase is presented. In figure- 
14, the TCA is powering-up from ignition to level-1 operation. As it did so, chug was exhibited with an 
amplitude of 10.3% of steady-state chamber pressure, with a spectral intensity of 14.424 psi/Hz 
(98.122 KPa/Hz) at a frequency of 195 Hz. The duration of this chug was about 139 milliseconds. 

In figure-15, the test-23 chamber pressure during the buzz-dominated phase is presented. In 
figure-15, the TCA was operating at level-2 conditions. As it did so, buzz dominated the spectrum 
with an intensity of 8.295 psi/Hz (56.429 KPa/Hz) at a frequency of 493 Hz. The buzz decayed away 
by about 35 seconds. In addition to the buzz, some low-intensity, low-frequency peaks in the 
spectrum appeared with a magnitude of 1.794 psi/Hz (12.204 KPa/Hz) at 59 Hz, 2.406 psi/Hz 
(16.367 KPa/Hz) at 127 Hz, 3.0 psi/Hz (20.408 KPa/Hz) at 195 Hz, and 1.505 psi/Hz (10.238 
KPa/Hz) at 269 Hz. The duration of the buzz during this phase was about 30 seconds. The maximum 
amplitude was about 7.6% of the steady-state chamber pressure. 

In figure-16, the test-23 chamber pressure during the buzz-minimized phase is presented. In 
figure-16, the TCA was still operating at level-2 conditions. As it did so, the intermediate-frequency 
oscillation appeared in the spectrum with a lower intensity of 1.147 psi/Hz (7.803 KPa/Hz) at a 
frequency of 498Hz. Additionally, some low-frequency oscillations appeared in the spectrum with a 
magnitude of 1.958 psi/Hz (13.320 KPa/Hz) at 59 Hz, 1.742 psi/Hz (11.850 KPa/Hz) at 181 Hz, 

0.490 psi/Hz (3.333 KPa/Hz) at 269 Hz, and 0.502 psi/Hz (3.415 KPa/Hz) at 298 Hz. The duration of 
this phase was about 114 seconds. The average amplitude was about 2.4% of the steady-state 
chamber pressure. 

The test-23 chamber pressure during the shut-down phase is presented. The TCA was 
powering-down from level-2 operation. As it did so, chug was exhibited in several sub-phases. The 
first and third sub-phases are presented in figures-17 and 18, respectively. In the first shut-down sub- 
phase, chug was exhibited with an amplitude of about 9.6% of steady-state chamber pressure, with a 
spectral intensity of 10.802 psi/Hz (73.483 KPa/Hz) at a frequency of 68 Hz, 3.831 psi/Hz (26.061 
KPa/Hz) at 132 Hz , and 3.828 psi/Hz (26.041 KPa/Hz) at 142 Hz. The duration of chug in the first 
shut-down sub-phase was 425 milliseconds. In the third shut-down sub-phase, chug was exhibited 
with an amplitude of about 18.4% of steady-state chamber pressure, with a spectral intensity of 
14.320 psi/Hz (97.415 KPa/Hz) at a frequency of 156 Hz. The duration of chug in the third shut-down 
sub-phase was 201 milliseconds. There was not a second shut-down sub-phase for test-23. 

Consider the test-24 chamber pressure during the start-up phase. Chug was exhibited with 
an amplitude of about 32.8% of steady-state chamber pressure, with a spectral intensity of 30.436 
psi/Hz (207.048 KPa/Hz) at a frequency of 195 Hz. The duration of this chug was about 185 
milliseconds. 

Consider the test-24 chamber pressure during the buzz-dominated phase. As the TCA was 
operating at level-2 conditions, buzz dominated the spectrum with an intensity of 13.635 psi/Hz 
(92.755 KPa/Hz) at a frequency of 498 Hz. The buzz decayed away by about 54 seconds. In 
addition to the buzz, some low-intensity, low-frequency peaks in the spectrum appeared with a 
magnitude of 2.259 psi/Hz (15.367 KPa/Hz) at 59 Hz, 5.983 psi/Hz (40.701 KPa/Hz) at 132 Hz, and 
3.211 psi/Hz (21.844 KPa/Hz) at 195 Hz. The duration of the buzz during this phase was about 48.5 
seconds. The maximum amplitude was about 13.3% of the steady-state chamber pressure. 

Consider the test-24 chamber pressure during the buzz-minimized phase. As the TCA was still 
operating at level-2 conditions, the intermediate-frequency oscillation appeared in the spectrum with 
a lower intensity of 2.130 psi/Hz (14.490 KPa/Hz) at a frequency of 503 Hz. Additionally, some low- 
frequency oscillations appeared in the spectrum with a magnitude of 1.749 psi/Hz (11.898 KPa/Hz) at 
59 Hz, 1.683 psi/Hz (11.449 KPa/Hz) at 181 Hz, 0.369 psi/Hz (2.510 KPa/Hz) at 273 Hz, and 0.512 
psi/Hz (3.483 KPa/Hz) at 298 Hz. The duration of this phase was about 98 seconds. The average 
amplitude was about 2.3% of the steady-state chamber pressure. 
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Figure-14: test-23 chamber pressure chug during the start-up phase. 
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Figure-15: test-23 chamber pressure during the buzz-dominated phase of mainstage. 
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Figure-16: test-23 chamber pressure during the buzz-minimized phase of mainstage. 
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Figure-17: test-23 chamber pressure during the shutdown first sub-phase. 
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Figure-18: test-23 chamber pressure during the shutdown third sub-phase. 


Consider the test-24 chamber pressure during the shut-down phase. As the TCA was 
powering-down from level-2 operation, chug was exhibited in several sub-phases. In the first shut- 
down sub-phase, chug was exhibited with an amplitude of about 12.4% of steady-state chamber 
pressure, with a spectral intensity of 13.809 psi/Hz (93.939 KPa/Hz) at a frequency of 78 Hz. The 
duration of chug in the first shut-down sub-phase was 404 milliseconds. In the second shut-down 
sub-phase, chug was exhibited with an amplitude of about 7.3% of steady-state chamber pressure, 
with a spectral intensity of 3.667 psi/Hz (24.946 KPa/Hz) at a frequency of 54 Hz and 6.203 psi/Hz 
(42.197 KPa/Hz) at 171 Hz. The duration of chug in the second shut-down sub-phase was 163 
milliseconds. In the third shut-down sub-phase, chug was exhibited with an amplitude of about 30.7% 
of steady-state chamber pressure, with a spectral intensity of 20.954 psi/Hz (142.544 KPa/Hz) at a 
frequency of 161 Hz. The duration of chug in the third shut-down sub-phase was 148 milliseconds. 

Consider the test-25 chamber pressure during the start-up phase. Chug was exhibited with 
an amplitude of about 29.5% of steady-state chamber pressure, with a spectral intensity of 31.203 
psi/Hz (212.265 KPa/Hz) at a frequency of 205 Hz. The duration of this chug was about 158 
milliseconds. 

Consider the test-25 chamber pressure during the buzz-free phase. There was no buzz- 
dominated phase. As the TCA was operating at level-2 conditions, some low-frequency oscillations 
appeared in the spectrum with a magnitude of 2.202 psi/Hz (14.980 KPa/Hz) at 59 Hz, 4.434 psi/Hz 
(30.163 KPa/Hz) at 142 Hz, 2.444 psi/Hz (16.626 KPa/Hz) at 181 Hz, 1.585 psi/Hz (10.782 KPa/Hz) 
at 215 Hz, and 1.522 psi/Hz (10.354 KPa/Hz) at 293 Hz. The duration of this phase was about 3.845 
seconds. The average amplitude was about 6.8% of the steady-state chamber pressure. 

Consider the test-25 chamber pressure during the shut-down phase. As the TCA was 
powering-down from level-2 operation, chug was exhibited in several sub-phases. In the first shut- 
down sub-phase, chug was exhibited with an amplitude of about 24.7% of steady-state chamber 
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pressure, with a spectral intensity of 37.995 psi/Hz (258.469 KPa/Hz) at a frequency of 78 Hz and 
9.805 psi/Hz (66.701 KPa/Hz) at 151 Hz. The duration of chug in the first shut-down sub-phase was 
263 milliseconds. In the second shut-down sub-phase, chug was exhibited with an amplitude of 
about 27.8% of steady-state chamber pressure, with a spectral intensity of 31.984 psi/Hz (217.578 
KPa/Hz) at a frequency of 195 Hz. The duration of chug in the second shut-down sub-phase was 
393 milliseconds. In the third shut-down sub-phase, chug was exhibited with an amplitude of about 
31.7% of steady-state chamber pressure, with a spectral intensity of 33.415 psi/Hz (227.313 KPa/Hz) 
at a frequency of 176 Hz. The duration of chug in the third shut-down sub-phase was 144 
milliseconds. 

Consider the test-26 chamber pressure during the start-up phase. Chug was exhibited with 
an amplitude of about 25.9% of steady-state chamber pressure, with a spectral intensity of 24.689 
psi/Hz (167.952 KPa/Hz) at a frequency of 186 Hz and 6.056 psi/Hz (41.197 KPa/Hz) at 376 Hz. 

The duration of this chug is about 187 milliseconds. 

Consider the test-26 chamber pressure during the first sub-phase of the buzz-free phase. As 
the TCA was operating at level-2 conditions, some low-frequency oscillations appeared in the 
spectrum with a magnitude of 2.424 psi/Hz (16.490 KPa/Hz) at 59 Hz, 0.903 psi/Hz (6.143 KPa/Hz) 
at 142 Hz, 2.417 psi/Hz (16.442 KPa/Hz) at 181 Hz, and 0.736 psi/Hz (5.007 KPa/Hz) at 298 Hz. The 
duration of this phase was about 48 seconds. The maximum amplitude was about 6.4% of the 
steady-state chamber pressure. 

Consider the test-26 chamber pressure during the second sub-phase of the buzz-free phase. 
As the TCA was still operating at level-2 conditions, some low-frequency oscillations appeared in the 
spectrum with a magnitude of 2.420 psi/Hz (16.463 KPa/Hz) at 59 Hz, 2.414 psi/Hz (16.422 KPa/Hz) 
at 181 Hz, and 0.431 psi/Hz (2.932 KPa/Hz) at 298 Hz. The duration of this phase was about 99 
seconds. The average amplitude was about 2.4% of the steady-state chamber pressure. 

Consider the test-26 chamber pressure during the shut-down phase. As the TCA was 
powering-down from level-2 operation, chug was exhibited in several sub-phases. In the first shut- 
down sub-phase, chug was exhibited with an amplitude of about 37.8% of steady-state chamber 
pressure, with a spectral intensity of 48.750 psi/Hz (331.633 KPa/Hz) at a frequency of 73 Hz, 

11.885 psi/Hz (80.850 KPa/Hz) at 151 Hz, and 5.445 psi/Hz (37.041 KPa/Hz) at 225 Hz. The 
duration of chug in the first shut-down sub-phase was 567 milliseconds. In the second shut-down 
sub-phase, chug was exhibited with an amplitude of about 9.9% of steady-state chamber pressure, 
with a spectral intensity of 24.486 psi/Hz (166.571 KPa/Hz) at a frequency of 63 Hz and 11.974 
psi/Hz at (81.456 KPa/Hz) 176 Hz. The duration of chug in the second shut-down sub-phase was 
126 milliseconds. In the third shut-down sub-phase, chug was exhibited with an amplitude of about 
36.2% of steady-state chamber pressure, with a spectral intensity of 27.069 psi/Hz (184.143 KPa/Hz) 
at a frequency of 156 Hz. The duration of chug in the third shut-down sub-phase was 159 
milliseconds. 


SUMMARY 

For the Fastrac TCA, tested at MSFC test stand 116, the injector pressure drops were 19% 
and 20% on the RP-1 and LOX sides, respectively. Using these pressure drops, a linear analysis was 
conducted to determine the low-frequency combustion stability characteristics of the TCA combined 
with the RP-1 and LOX feed systems downstream of the venturies. This analysis was only valid when 
making predictions of low-frequency combustion stability during mainstage operation only. The most 
obvious result of this analysis was that the chamber pressure oscillations were stable. The low- 
frequency combustion stability was not very sensitive to changes in AP F /Pc- According to figure-8, 
APf/Pc could be lowered to below 4% without chug being exhibited. However, The low-frequency 
combustion stability was more sensitive to changes in AP 0 /Pc. According to figure-8, AP 0 /Pc could be 
lowered to about 10% without chug being exhibited. Unfortunately, special testing to verify the 
sensitivity of low-frequency combustion stability to injector pressure drops was never performed. 
Additionally, these pressure oscillations would occur in the TCA and RP-1 and LOX feed systems at a 
frequency of 65 Hz and would be damped to 50% of their initial amplitude in about 10 milliseconds. 
The maximum allowable damp time was about 155 milliseconds. 
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For the Fastrac flight engine, tested at SSC, the injector pressure drops were unchanged at 
19% and 20% on the RP-1 and LOX sides, respectively. Using these pressure drops, the linear 
analysis was repeated to determine the low-frequency combustion stability characteristics of the flight 
engine consisting of the TCA combined with the RP-1 and LOX feed systems. For the flight engine, 
the propellant feedlines were shortened according to flight engine design in figures 9 and 10. 

Additionally, in the analysis, the venturis were not replaced by pumps. The result of this analysis was 
that the chamber pressure oscillations were stable. According to figure-1 1 , the low-frequency 
combustion stability was insensitive to changes in APf/Pc above AP 0 /Pc = 10%. Also, low-frequency 
combustion stability was insensitive to changes in APo/Pc below APf/Pc = 20%. Unfortunately, the 
frequency and damp time of these pressure oscillations that would occur in the TCA and in the RP-1 
and LOX feed systems of the flight engine were indeterminable. 

In order to convert from a pressure-fed TCA to a pump-fed TCA, at the least, the propellant 
feedlines were shortened. The effect of shortening the propellant feedlines on low-frequency combustion 
instability was illustrated in figure-12. Shortening the propellant feedlines caused the neutral stability 
boundary to curve downward for low values of APf/Pc, thereby increasing the region of stable 
combustion. In the following appendix, equations 13, 33, and 36 indicates that decreasing the propellant 
feedline volumes results in “stiffening” the propellant feed system. This is a stabilizing effect. 

The experimental data from tests 23 through 26 indicated that during mainstage, low-frequency 
chamber pressure oscillations were present. Typical frequencies consisted of the fundamental frequency 
of 60 Hz with its first and second harmonics at 180 Hz and 300 Hz, respectively. An additional group of 
peaks in the spectrum occurred at frequencies of 127 Hz, 132 Hz, and 142 Hz. A second group of 
peaks in the spectrum occurred at the frequencies of 195 Hz and 215 Hz. Finally, a third group of peaks 
in the spectrum occurred at the frequencies of 269 Hz and 273 Hz. However, the spectral intensity of 
these peaks were greatly overshadowed by the peak that was at 500 Hz in tests 23 and 24. Once the 
500 Hz intermediate frequency combustion instability was eliminated in tests 25 and 26, the low- 
frequency pressure oscillations remained but resulted in low amplitudes that did not raise any concerns 
about combustion instability during mainstage. 

The experimental data from tests 23 through 26 indicated that during start-up, chug was present. 
The duration of the chug was from 139 to 187 milliseconds. In the chamber pressure frequency 
spectrum, there was a peak at a frequency of 186 Hz. Also, there was a typical group of peaks which 
had frequencies at 195 Hz and 205 Hz. Finally, there was a peak with a frequency of 376 Hz. Note that 
the peak at 376 Hz was secondary in spectral intensity and was very close to the first harmonic of 186 
Hz. 


The experimental data from tests 23 through 26 indicated that during shutdown, chug was 
present. The duration of the chug was from 601 to 852 milliseconds. In the chamber pressure frequency 
spectrum, there were peaks at the frequencies of 78 Hz, 151 Hz, and 225 Hz during the first shutdown 
sub-phase. Note that the peaks at 151 Hz and 225 Hz are close to the first and second harmonics of 78 
Hz, respectively. Also, during the first shutdown sub-phase, there were peaks at the frequencies of 68 
Hz, 132 Hz, and 142 Hz. Note that 132 Hz is close to the first harmonic of 68 Hz. During the second 
shutdown sub-phase, there were a group of peaks at the frequencies of 54 Hz and 63 Hz. Note that 
these frequencies are close to 60 Hz. Also, during the second shutdown sub-phase there were another 
group of peaks at 171 Hz and 176 Hz. Note that these frequencies are close to 180 Hz, the first 
harmonic of 60 Hz. Finally, during the second shutdown sub-phase there was a peak at 195 Hz. During 
the third shutdown sub-phase there were peaks at the frequencies of 156 Hz, 161 Hz, and 176 Hz. 

While the analysis indicated damped pressure oscillations at 65 Hz, the test data indicated peaks 
at low spectral intensity at 60 Hz. This leads to one of two possibilities. One possibility is that the analysis 
accurately predicted a frequency of pressure oscillation. The other possibility is that this close agreement 
is merely coincidental. The appearance of the 60 Hz oscillations with the respective harmonics is due to 
poor transducer cable connections resulting in “ground looping 7 ” or poor grounding of the cable. The 
causes of ground looping may be numerous. One cause may be moisture trapped between the 
transducer and cable lead. 


19 


ACKNOWLEDGMENTS 


The authors would like to acknowledge and thank Dr. John J. Hutt (NASA/MSFC), team leader of 
the Fastrac Combustion Devices Team, whose contributions were crucial during the design of the Fastrac 
injector and whose direction was crucial during the Fastrac TCA performance testing. Dr. Hutt originally 
requested the development of this analysis. 

REFERENCES 

1) “Guidelines For Combustion Stability Specifications And Verification Procedures For Liquid 
Propellant Rocket Engines”, CPIA Publication-247, page-4, 1973. 

2) “Elimination of High-Frequency Combustion Instability in the Fastrac Engine Thrust Chamber", 

M. Rocker and T.E. Nesman, presented at the 1999 JANNAF Joint Meeting of the Combustion, 
Airbreathing Propulsion, and Propulsion System Hazards Subcommittees in Cocoa Beach, Florida, 

1999. 

3) “Elimination of Intermediate-Frequency Combustion Stability of the Fastrac Engine Thrust 
Chamber”, M. Rocker and T. E. Nesman, to be published, 2001. 

4) “Analysis of Chugging in Liquid-Bipropellant Rocket Engines Using Propellants with Different 
Vaporization Rates”, L.M. Wenzel and J.R. Szuch, NASA TN D-3080, 1965. 

5) “Rocket Propulsion Elements”, fourth edition, G.P. Sutton and D.M. Ross, pages 276-277, 

1976. 

6) “Liquid Propellant Rocket Combustion Instability”, NASA SP-194, editors: D.T. Harrje and F.H. 
Reardon, page-269, 1972. 

7) Private Communication, R.A. Hewitt, Aerojet, June 14, 2001. 

APPENDIX: REVIEW OF MODEL FORMULATION 

The formulation of the model by which low-frequency combustion instability predictions were 
performed began with the modeling of the liquid propellant feed system. A feed system element was 
defined, as illustrated in figure-19. The feed system element had two components, a feedline 
segment and a flow control device. The feedline segment was represented by a volume of fluid. 

While the pressure within the feedline segment was assumed to vary with time, it did not vary in 
space. Additionally, the feedline segment received fluid from upstream and lost fluid through the flow 
control device downstream. The flow control device was represented by a flow constriction between 
two feedline segments. It was understood that the flow control device could have been an orifice, 
valve, or anything across which a pressure drop occured. Within the flow control device, the flow rate 
varied with time but was constant along the length of the flow control device. 


d 



feedline j flow | 
control 
device 


Figure-19: feedline element schematic. 


20 



_d l_ 
dt 


Within the feedline segment, the fluid behavior was governed by the conservation of mass: 

IjifpdV + |j>pv»Wds = 0 . 


( 5 ) 


After applying conservation of mass to the feedline segment in figure-19, equation-5 became: 


V T 


dPi _ 


dt 


= m a - m 


Fluid density is usually presented as a function of pressure and temperature: 

A = p(p i-Ti) • 

However, since isothermal flow was assumed, equation-7 became: 

A = P( Pi) • 

Therefore, the term on the left-hand side of equation-6 became: 
dp, _ dp, dP, _ ^dE> 


( 6 ) 


(7) 


( 8 ) 


dt 


dP, dt 


p'( p .) 


dt 


(9) 


Since the fluid was assumed to be a liquid at constant temperature, density was assumed to vary 
linearly as a function of pressure. Therefore: 


p'( ) = constant 


( 10 ) 


Therefore, substituting equation-9 into equation-6 and linearizing about steady-state conditions 
gave: 


V LINE P'( Pj 


dt 


= m, - m 


( 11 ) 


Scaling equation-11 by the steady-state flow rate and steady-state chamber pressure gave: 

( p ,N l 


cA 

dt 


V P c7 
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u 

m 


m 

m 


( 12 ) 


Where the scaled capacitance was defined as: 
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c = 


^VlinePIPi) 


m 


( 13 ) 


It was determined that the scaled capacitance had the units of seconds. Note that for the steady- 
state flow rates: 


= m 


(14) 


Within the flow control device, fluid behavior was governed by the conservation of momentum 
of Newton's second law of motion: 


dt 


|£fpVdV + fypv(v • N)ds = -^PWds + Jj\ds . 


Where for the first term on the left-hand side of equation-15: 
d rrr d r „ „ dm 


-§pvd v = -JpuAdx e x = 


dt- 


dt 


For the second term on the left-hand side of equation-15: 


§pv{VN ) ds = (p e u 2 e -p iU ^)A e x = 


1_ 

VPe 


Pi 


m 2 

A 


For the first term on the right-hand side of equation-15: 

£j>PWds = (p e - pJa e x . 

And, for the second term on the right-hand side of equation-15: 


IK 


ds = 


f PavgKvg ^ & 


f 1 f m 


4 2p av AA 


Ind e. 


(15) 


(16) 


(17) 


(18) 


(19) 


By substituting equations-16, 17, 18, and 19 into equation-15, conservation of momentum has been 
applied to the flow control device in figure-19. Therefore, equation-15 became: 


dt 


r 1_ 
Ae 


Pi 


m 


= (p ± - P e )A - I tnd 

A Vl 4 2p V. A, 


( 20 ) 


The cross-sectional area of the flow control device and the average density across the flow control 
device was defined, respectively, by: 
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A 


7Td 2 

— ' 2 Pavg = Pi + Pe • 


(21), (22) 


Dividing equation-20 by the cross-sectional area, and then substituting in 


equations 21 and 22 gave: 


£ dm 
A dt 



IPs 



££ 1 f m 

d (p. + p e )U 


The flow control device loss coefficient was defined such that: 


(23) 


K 


loss 


— : fanno line 
d 

Kioss • orifice, valve, etc. 


Therefore, equation-24 became: 


£ dm 
A dt 


r 1 

vPe 


Pi 


— = P, 


Kioss ( m 


(Pi + Pe)v A 


(24) 


(25) 


Linearizing equation-25 about the steady-state conditions gave the ordinary differential equation for 
the perturbation in flow rate: 


£ dm' 
A dt 


1 - 



1 - 


p'( p .) 


pi 


^ rh ^ 




K 


loss 


] 2mm' 


(Pi + Pe)J A 2 


K, 


.P'( p i) 


(Pi + Pe) 

K 


rmV 


p' 


loss- 


pp.) 

2 


(Pi + Pe) 


m 

VAy 


p: 


(26) 


In the process of linearizing equation-25, the steady-state momentum equation was obtained: 


K 


loss 


m 


4 


(Pi + Pe) V A, 


= Pi - Pe - 


1 

iPe 


Pi) 


f —\ 2 

m 

vAy 


(27) 


Substituting equation-27 into equation-26 and dividing through by the steady-state chamber 
pressure gave: 


L 


— f— 1 

dtV m ) 



m 



( 28 ) 
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Where the scaled inertance and the scaled resistance were respectively defined as: 


L = A! _ R = 2 ( P , _■ P .) 

P C A ' P c 

(29), (30) 


The units for the scaled inertance and scaled resistance were seconds, and dimensionless, 
respectively. The inlet convective-resistance was defined as: 


K; 


= 1 - 


P^i) 


Pi 






pM 


(Pi + Pe) 


P - P. - 


Also, the exit convective-resistance was defined as: 


vPe 


_1 

Pi) 


V m^ 2 


vAy 


(31) 


K = 1 


P'( P e) 




Pe Uy 


Jp 

(ft * ft) ‘ 


R - P, - 


1 

Pe 


Pi 


i J 


Ay 


The units for both the inlet and exit convective-resistance were dimensionless. Differentiating 
equation-28, and substituting in equation-12 gave: 


(32) 


d 2 


no 

_ d 



K, m' 

K. m' 

d 

f P e '"j 


+ R — 


+ 


l_ U 

- K 

e_ 

V m ) 

1 dt v 

v m J 


C rh 

C m 

dt 

^ J 


dt 

Combined with equation-12: 


dt 


S' 

v p cy 


rh^ 

m 


m' 


m 


(33) 


(34) 


these two ordinary differential equations describe the fluid flow inside of the feed system element. 


Equation-33 is dynamically similar to two analogous systems. The first system is electrical and 
the second system is mechanical. For the electrical system, the governing equation was: 

_ d 2 Q „ dQ 1 dE 

L — + R— + -Q = — 
dt 2 dt C dt 

(35) 


Where Q was the electronic flow rate or current. By comparing equation-35 to equation-33, one 
concluded that L was either the inductance in equation-35 or the scaled inertance in equation-33. 
Also, R was the scaled resistance in equation-33 and the resistance in equation-35. Finally, C was 
the scaled capacitance in both equation-33 and the capacitance in equation-35. For the second 
system, the governing equation was: 

d 2 x dx , . x 

m — ^ + c — + kx = f(t) 

dt 2 dt v ' 

(36) 
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Where x was the displacement of mass-m. By comparing equation-36 to equation-33, one concluded 
that m, which was the linear inertia or mass in equation-36, was analogous to the inertance-L in 
equation-33. Also, c, which was the damping coefficient in equation-36, was analogous to the scaled 
resistance-R in equation-33. Finally, k, which was the stiffness in equation-36, was analogous to 1/C, 
where C was the scaled capacitance in equation-33. A large scaled resistance in equation-33 would 
contribute to the elimination of oscillations in the feed system. Since the stiffness in equation-36 is 
inversely proportional to the scaled capacitance in equation-33, A feed system with a large scaled 
capacitance will have a low stiffness. Such a feed system was said to be “soft”. A soft feed system 
will have either large feedline volumes, feedline volumes filled with highly compressible fluids, or both. 
A feed system with a low scaled capacitance will have a high stiffness. Such a feed system is said to 
be “stiff’. A stiff feed system will have small feedline volumes, be filled with a low compressibility fluid, 
or both. 


rh'(t) 

P'(t) 


It was assumed that the solutions of equations-33 and 34 would be exponential functions: 
= M(s)e st , m;(t) = M„(s)e st , 


= P i (s)e st , P e '(t) = P e ( s)e 


st 


(37), (38) 


(39), (40) 


The exponential function exhibited sinusoidal behavior with damping when s was defined as: 

s = X + jco , j = V-l . 


(41), (42) 


Substituting equations-37, 38, 39, and 40 into equations-33 and 34 gave: 


kAm{s) 


LS + RS + 

m 


Kj M a (s) 
C rh 



Cs 


PM 

Pc 


MjM _ M(s) 

m m 


(43) 


(44) 


After the equations of fluid motion had been established for a single feed system element, 
consisting of a feedline segment and a flow control device, the equation set describing the dynamics 
of a single propellant feed system consisting of a series of several feed system elements was 
established. The model of a single propellant feed system is presented in figure-20: 



combustion 

chamber 


Figure-20: single propellant feed system schematic. 


25 




The three feedline elements linked in series model the portion of the feed system 
downstream of the venturi and upstream of the combustion chamber. This resulted in a set of three 
momentum equations, one for each of the three elements: 


( 


LjS + R x s + 


L 2 s + R 2 s + 


L3S + R3S + 


K; 




'1 J 


K 


i2 


c 


Mi(s) = K n M v _(s) 
m C 1 m 

M 2 (s) _ K i2 M 2 (s) 


K . S 


P 2 (s) 


el 1 


2 J 


m 


C 2 m 


K e2 s 


^(s) 


P c 


K i3 \M 3 (s) 


_ K i3 M 2 (s) 


'3 J 


m 


C 3 m 


K e3 s 


P c(s) 


P c 


Where for the mass equations: 


C 2 S — 


P 2 ( s) _ ^(s) 


M 2 (s) 


m 


m 


C 3 s 


-Pj(s) _ M 2 (s) M 3 ( s) 


m 


m 


Equations 48 and 49 were substituted into equations 45 and 46, respectively. The equations 
became: 


a,i(s) a n 0 


'^(s) 


0 

a 2\ a 22( S ) a 23 


Y 2 ( S ) 

► = 

0 

0 a 32 « 33 (s) 


Y 3 (s) 


yM. 


Where for the cc’s: 

°ai( s ) = 


+ R 3 s + 


K. 


K, 


C, 


(45) 


(46) 


(47) 


(48) 


(49) 


(50) 


(51) 


a. 


K 


el 




K 


i2 


(52), (53) 


^22(^0 L 2 s + R 2 s + 


K, 


K. 


(54) 
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(55), (56) 


a 


23 



a 


32 



a 33 (s) = 


L 3 s + 


R 3 s + 


K 


i3 


c. 


Also, for the admittances: 


Yi(s) 


Mi Ml 

P c (s ) m 


Y 2 ( S ) 


M 2 (s) P c 

p c ( s ) m 


Y 3 (s) 


^( s ) § 

P c (s) m 


(57) 


(58), (59), (60) 


Finally, for the chamber pressure perturbation forcing function: 

X 3 ( S ) = " K e3 S • 


Note that for the flow rate perturbation through the venturi: 

M v ( s) = 0 . 


Solving equation-50 for the admittances gave: 


Y 3 (s) = 


a 


^11^23^32 


33 


(^21^12 ^11^22) 


(61) 


(62) 


(63) 


The ex’s and y 3 were assumed to correspond to the oxizider feed system. Equation-63 became: 


Y 30 (s) 


O' 


r 3 

^ 11 ^ 23^32 

{ a 2i a n ~ ^11^22) 


Y 0 (s) • 


(64) 


Similarly, the (3’s and 63 were assumed to correspond to the fuel feed system. Therefore, equation-65 
gave: 


Y 3F (s) 


A 


A1A1A 


32 


(A1A2 - P1A2) 


y p (s) 


(65) 


In figure-10, schematic of the combustion chamber is presented. First, the propellants 
undergo injection into the chamber. Second, the propellants are atomized, vaporized, mixed, and 
undergo combustion or burning to produce combustion products. Third, the combustion products are 
expanded through the nozzle. In the combustion chamber, the combustion products were assumed 
to behave as an ideal gas. Additionally, the pressure and temperature in the combustion chamber 
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behaved as a function of time only, and did not vary in space. Therefore, conservation of mass with 
the nozzle gas flow rate gave: 


Vc 

R 


_d_ 

dt 


( 





<(t) 


m N (t), m N 


1 y + 1 


2U-1 


Y + 1 



(66), (67) 



injection | combustion . expansion 
or burning 


Figure-21: schematic of rocket engine combustion chamber. 


Where for the propellant burning rate and the mixture ratio: 


iru 


^FR ' 


rru 


(68), (69) 


Substituting equations 67, 68, and 69 into equation-66 and the linearizing equation-66 about the 
steady-state conditions gave the ordinary differential equation for the perturbations in chamber 
pressure and chamber temperature: 


P c v c ^ 

fid 

Pc V c d I 

X C 

m dt 

UJ 

m dtl 

UJ 


f r , 

f i y™ 

f p' 
_c_ 

1 T c ' 

Vr + lj m 0B 

Vr + lj m FB 

V^C 

2^ 


\ 

J 


(70) 


The chamber temperature is a function of both chamber pressure and mixture ratio. However, it is 
regarded as being a strong function of mixture ratio only: 


T c = T c (r) 


(71) 


The Taylor’s series expansion of the chamber temperature about steady-state conditions gave: 
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( 72 ) 


T c (r) = T c (r + r') = 


\ ,<?r r . 1 / A 2 0 » 2 T r 

T ' (r) + + k r) + 


Linearizing gave: 


T c (r) = TJ. 


+ T' = 


/— \ / <?Tr 

T c( r ) + 


Therefore, for the perturbation in chamber temperature: 


r £?T c r' 


dr r 


From the linearization of equation-69, perturbation in mixture ratio was determined to be: 


(73) 


(74) 


r 

r 




rru 


m c 


itu 


(75) 


Substituting equation-75 into equation-74 gave: 


_c _ 

% 


for 


r ■ / 

m r 


^ <3r 


m c 


V ^OB 


m 


FB J 


(76) 


Substituting equation-76 into equation-70 gave: 


p c v c a 

f pO 

fc_ 

, P C _ PcV c r dT c d 

( ■' \ 
m OB 

p c V c r r?T c d 

f m' ^1 

iU FB 

m dt 


P c m TJ. dr dt 

\ ^OB ) 

m TJ, dr dtJ 

v ^ fb y 


fr * ) 

| + 1 r dr c ] 

, ™OB I 

1 1 

| Irk' 

[v r + lj 

1 ‘ 2^ dr 

k, i 

\r + lj 

2% dr 


rru 


nu 


(77) 


For the perturbations in the oxidizer and fuel burning rates: 

™OB(t) = ' Kb^) = m;(t-T p ) . 


(78), (79) 


Equations 78 and 79 stated that the burning rate of a given propellant at some instant t was not 
equal to the injection flow rate of the given propellant at the same instant. But the burning rate at 
some instant t was equal to the injection flow rate at the prior instant t-io or t-x F for the oxidizer or 
fuel, respectively. For the derivative with respect to time of the perturbations in the propellant flow 
rates: 
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(80), (81) 


dt 

dt 


d lh °( t " T °)^( t_ ‘ r o) 


d (t-T 0 ) 

d 


— 7 rltl'ft - T f ) (t - T„) 

d(t-T F ) FV F 'dt v F) 


d 


d (t-r 0 ) 

d 

d(t-T F ) 


<(t-r 0 ) , 
K{ t-Tp) ■ 


For the combustion chamber flush/fill time, the steady-state oxidizer flow, fuel flow, and total 
propellant rates: 


Me 

m 


, m r 




m F , m — m n + m F 


(82), (83), (84), (85) 


Substituting equations 78 through 84 into equation-77 gave: 


d 

T 


P' r <?T r 

1 L — t e 

d 

*0 

t-gi 

* C 

dt 

Uj 

‘ — * r — ^ 

P c ^ dr d 

r <?T C 

T u 

O 

1 

l 

<( 

*0 j 


c TJ, dr d( 

n r 
1 

^ J 


r + 1 
1 

r + 1 


1 r ^T c ] m 0 (t-T 0 ) 

2 % dr J m 0 

1 r m'(t - r F ) 

2 dr m F 


Where for the perturbation in chamber pressure: 


P'(t) = P c ( s)e 


st d , , st 


dt 


p c(t) = sP c( s K 


For the perturbation oxidizer flow rate: 
A'o(t) = M 0 ( s)e St , 

(t-To) 


= M 0 (s)e 


' d(t - T 0 ) 


™o(t- T 0 ) = sM 0 (s)e 


And, for the perturbation fuel flow rate: 

m F (t) = M F (s)e st , 

' F (t-r F ) = M f ( s)e 


m 


s t-r. 


d 


' d(t - r F ) 


m F (t-r F ) = sM F (s)e 


( 86 ) 


(87), (88) 


t t 0 j 


(89), (90), (91) 


s[t-T F ) 


(92), (93), (94) 
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Substituting equations 87, 88, 90, 91, 93, and 94 into equation-86 gave the characteristic equation: 


r c s + 1 = f 0 (s)Y 0 (s)e ST ° + fp(s)Y F (s)e STr . 

(95) 


Where the injector admittances for the oxidizer and fuel were given by: 

, x mJs) P r , s mJs) P r 

Y 0 (s) = , Y f (s) = fW c 


P c (s) rh c 


For the oxidizer and fuel gains: 


PM 


m. 


(96), (97) 


f o( s ) = 


f F (s) = 


r + 1 
r 

r + 1 
1 

r + 1 
( 1 


: constant gas temperature assumed 


+ 1 T C S + — 


13 r <?T„ 


: fluctuating gas temperature assumed 

(98) 


2)\dr 

: constant gas temperature assumed 


r + 1 


-I T c s + - 


13 r 


2) % dr 


: fluctuating gas temperature assumed 


(99) 


Dividing equation-95 by the left-side and then multiplying equation-95 by the denominators of both 
Yo(s) and Y 0 (s) as given by equations 64 and 65, gave: 


(ctjj + A )(/?33 + B j f 0 C E (/? 33 + B j + f F D F (<Y 33 + A ) 


(100) 


Where: 

A* = 


Cl (I Cl 

^ 11 ^ 23^32 


^ 21^12 ^ 11^22 


, B = 




P21P12 - P11P22 


( 101 ), ( 102 ) 


-stt 


■ST, 


C* = / 3 , D* = <5 3 , E* = 


r c s + 1 


, F = 


T C S + 1 


The characteristic equation, equation-100 was re-written as: 


^33^33 + aa 33 + ^ 


33 


c . 


(103), (104), (105), (106) 


(107) 


Where for the coefficients a, b, and c: 
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a = B* - f P D*F* , b = A* - f 0 C*E* , 
C = f 0 C*B*E* + f p D*A*F* - A*B* . 


(108), (109) 
( 110 ) 


The characteristic equation, equation-107 was split into the real and imaginary portions. For the real 
portion: 


^■1^30^3F ^-'2^30 ^3^3F ^4 ® 


For the imaginary portion: 

^l^O^F + b 2^30 + b 3^3F + F> 4 — 0 . 

Where for the coefficients Ci, C 2 , C 3 , and C 4 : 

C P = A 2 - ry 2 , 

C 2 = (A^ F -ror 7 F ) + (a R A-a r <») , 

c 3 = (^o- W7 ?o) + ( b R A " b I W ) ' 

^“4 — ('SO^F — VoVf) + (^R^O - ^iVo) + ( b R^F — D I^7 f) - C R 

And for the coefficients Di, D 2 , D 3 , and D 4 : 

D P = 2 Xcl) , 

D 2 = (At7 F +ft>4) + (a R a> + a P A) , 

°3 = (°^o +Ar lo) + ( b R W + b I A ) - 

b 4 — ('SO^F + Vo^f) + {^rVo + ) + ( b R^?F + b I*3F ) ~ *"1 

Where from equations 108, 109, and 110: 

a = a R + ja 3 , b = b R + jb P , c = c R + jCj. . 

For the oxidizer feed system: 


(111) 

( 112 ) 

(113) 

(114) 

(115) 

(116) 

(117) 

(118) 

(119) 

( 120 ) 

(121), (122), (123) 


£ 0 = L 3o(^ 2 -<y 2 ) + ^ < % = 2l 30 aw . 

b 30 


(124), (125) 


For the fuel feed system: 
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(126), (127) 


£ f = L 3? (A 2 -(0 2 ) + , % = 2L 3F Aft) 

C 3f 


Solving equations 111 and 112 for R 3 o gave: 


R 
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c 3 r 3f 

+ o 


^3^3F 

+ d 4 ] 

V ^1^3F 

+ C 2.J 


V ^ > 1^3F 

+ D 2y 


Solving equation-128 for R 3F gave: 

EiR 3f + ^2^3F + ^3 — 0 / 

("E 2 ± a/ E 2- 4E : E 3 ) • 


R 


1 


2E 1 

Where for Ei, E 2 , and E 3 : 



= c 3 d 3 

- CjD- 

E 2 

= c 3 d 2 

- c 2 d. 

E 3 

- c 4 d 2 

- C 2 D, 


CiD 4 + c 4 d x 


(128) 


(129) 

(130) 

(131) 

(132) 

(133) 


Once R 3F was determined, it was substituted into equation-128 to determine R 3 o- The scaled 
resistances R 3 o and R 3F were defined as: 


R, 




(134), (135) 

The scaled resistances given by equations 134 and 135 are twice the scaled pressure drops across 
the oxidizer and fuel injectors elements. For a given frequency range and single value of damping 
coefficient, a curve of R 30 /2 vs. R 3F /2 may be generated. For a non-zero damping coefficient, if the 
curve passes through the operating point of the injector, then the frequency will correspond to the 
chamber pressure oscillations and the damping coefficient will give the damp time. A negative 
damping coefficient will correspond to chamber pressure oscillations that will decay with time. A 
positive damping coefficient will correspond to chamber pressure oscillations that will grow with time. A 
zero damping coefficient will correspond to the neutral stability boundary, where chamber pressure 
neither grow or decay. 
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